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. Abstract 

S ■ Movements of molecular motors on cytoskeletal filaments are described by directed 

, walks on a line. Detachment from this line is allowed to occur with a small prob- 

^ I ability. Motion in the surrounding fluid is described by symmetric random walks. 

Effects of detachment and reattachment are calculated by an analytical solution of 
the master equation. Results are obtained for the fraction of bound motors, their 

I average velocity and displacement. Enclosing the system in a finite geometry (tube, 

■ slab) leads to an experimentally realizable problem, that is studied in a continuum 

. description and also numerically in a lattice simulation. 
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I ■ 1 Introduction 
^ ■ 

c : 
o . 

O ' The interior of cells is both highly structured and dynamical. Active transport 

^ ■ therefore plays a crucial role to target molecules to the various compartments 

^ . of the cell as well as to maintain and reorganize the cell's structure. These 
tasks are accomplished by the cytoskeleton, a network of protein fibers, and 
the cytoskeletal motors, proteins which use the filaments of the cytoskeleton 
as highways or rails for directed transport. 



O 
O 



Cytoskeletal motors are enzymes which catalyze the hydrolysis of adenosine- 
triphosphate (ATP) and, at the same time, use the free energy released from 
this reaction to perform mechanical work and convert it into directed move- 
ments. These motor molecules have been studied extensively during the last 
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decade, both experimentally and theoretically. The main emphasis has been 
on the properties of single motor molecules, their motor mechanisms, and their 
directed walks along filaments [1,2]. 

The nanometer size and piconewton forces of the molecular motors imply that 
the typical binding energies are of the order of the thermal energy and 
that the motor-filament binding can be overcome by thermal fluctuations. On 
large time and length scales Is and 3> l/xni), motors perform peculiar 
random walks, where periods of directed active movements along filaments 
alternate with periods of non-directed Brownian motion in the surrounding 
fiuid after unbinding from a filament [3,4,5,6]. (A related problem is the effec- 
tive diffusion of an adsorbed particle along a surface via bulk excursions [7].) 
In order to study these random walks, we have recently introduced a class of 
lattice models, which are, on the one hand generic in the sense that they are 
independent of the specific motor mechanisms, but which we can, on the other 
hand, also apply to describe specific motor molecules by adapting the model 
parameters to the observed transport properties [4]. In addition, we can eas- 
ily incorporate motor-motor interactions such as the mutual exclusion from 
binding sites of the filaments, which leads to 'traffic jam'-like density pattern 
and various kinds of phase transitions [4,8,9] 

In the following, we will review our results for the motor's random walks 
obtained from these lattice models and present some new results from the 
corresponding continuum equations. The article is organized as follows: We 
introduce the lattice models in section 2. In section 3, we summarize our 
analytical results for the random walks on a lattice without confining walls. 
Compartments with confining walls which are experimentally accessible are 
studied in section 4, where we summarize our numerical results and present 
new analytical results from a continuum description of the random walks. 



2 Lattice models 

We have studied the random walks arising from many encounters with fila- 
ments by mapping them to random walks on a lattice [4,5]. A line of lattice 
sites represents a filament. Motors at these sites perform a biased random walk 
and move predominantly in one direction. Per unit time, they make a forward 
step, a backward step and no step with probability 1 — 7 — 6/2 — e{d — l)/d, 
5/2, and 7, respectively, where d = 2, 3 is the spatial dimension. With a small 
probability e/2d, they move to each of the adjacent non-filament sites and thus 
unbind from the filament. At the non-filament sites the motors perform simple 
symmetric random walks and move to each neighbor site with probability l/2d 
{d denotes the spatial dimension) and rebind to the filament with probability 
TTad when they reach again a filament site. Confining walls are implemented as 
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repulsive boundaries, at which all movements into the walls are rejected. 



3 The motors' random walks on a lattice 



Consider a discrete time random walk on a two dimensional square lattice with 
lattice sites labeled by integer coordinates (n, m) with the above transition 
probabilities. The master equation for this dynamics reads [5,6] 



Pn,m{t ~l~ 1) ^ Pn+l,m ~l~ ^ Pn—l,m ^ Pn,m+1 ~l~ ^ Pn,ra—1 {jTl ^ 0, il) 

Pnfi{t + 1) = ^ + \ Pn,-1 + (1 - 7 " ^)Pn-l,0 + ^Pn+1,0 + ^Pnfi 

Pn,±l{'t + -'-) ^ ^ -Pn+l,±l + ^ -Pn-l,±l + ^ Pn,±2 + '^Pnfl- (1) 

As initial condition we take an ensemble of particles at n = m = 0. Note that, 
for simplicity, we have chosen the sticking probability TTad = 1- 

The Fourier-Laplace transforms of the probability distribution along the fila- 
ment Ph{n,t) = Pn,o{t) and of the full distribution Pn,m{t) are defined as 



CXD CXD 



'irn p ( j-\ CXD CXD ■Kjiii + i rri p 

t=0 n=-oo \^ ^ ''J t=0 m,n=-oo \^ ^ ''J 



The master equations is reduced to an algebraic equation 



P{q,r,s) 



1 + [7(1 — cosr) + i2^(cosr — cosg) + iv\) sin r]Pb{r, s) 



s + 1 — I cos q 



cosr 



(3) 



where v\^ = l — ^ — 6— ^eis the average speed on the line. By integrating this 
result over q we also obtain Pb(^5 -s) on the left hand side. It thus satisfies a 
linear equation, that is easily solved. We end up with the distribution 



Pb{r,s) — [s + {l — '-f){l—cosr) + ^e{cosr — e — it'b sin r] ''■ (4) 

for the motors bound to the filament, where cosh/x = 2 + 2s — cosr. The 
distribution for the unbound motors follows from Eq. (3). 
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3. 1 Properties of the motors hound to the filament line 

Survival fraction. We extract the transport properties of the motor's random 
walk from (4). The value at r = gives the Laplace transform Nq{s) of the 
probability A^o(i) = Sn-Pn,o(i) that the motor is bound to the filament: 



For small times this implies A^o(^) ~ 1 — 2 e-\/I/ y^+e^t, which is a power series 
in e\/t. For t <^ this is somewhat surprising: although the motors detach 
at times ~ 1/e, the recurrent behavior of the random walk brings them mostly 
back to the filament, with half-integer powers in t due to diffusion. For times 

1 > l/e^ we find No{t) [1 - l/{2eH)]/\/'K e^t. The ^jg^ay expresses 
that finally all motors unbind, and agrees with scaling arguments [3,4]. 

Average position and speed on the filament line. The expression (4) for P\,{r\ s) 
contains much more information. At linear order in r we may derive the av- 
erage position of motor particles along the filament fine, Ni{t) = J2n'>^Pn,o{'t)- 
We obtain for short times A^i(t) = v^t [1 - 8e^/i/{3^/n)]. The average position 
of the motors bound to the filament is given by nb(^) = Ni(t) / NQ(t) ^ Vi,t{l — 

2 e\/t/3\/7r) and their average speed is v\y = dfib/dt ~ fb(l — e\/i/ \fn), where 

is the average speed if the particles did not leave the line. For large times one 
gets nio{t) ■UbV^fl — 2/(e\/7rt)]/e and %{t) ^ V\,^l1e\/t — (7r/2)A^o(*) '^^b, 
confirming the scaling v\^{f) ~ fb-^o(^) [3,4]. The effective motor velocity is 
reduced by a factor ~ -^o(t), i.e., by the probability that a motor is in the 
bound state. The relation ~ A^o^^b also applies to a simple two-state random 
walk, where motion is directed in one of the states only. In contrast to the 
simple two-state random walk, however, the probability A^o is time-dependent. 

Dispersion and diffusion coefficient on the filament line. From the second 
moment A^2(^) = Z)n ^^-Pn,o(^) we may define the normalized second moment 
n^b = -^2(^)7 Ao(^), the dispersion An^ = n^^y—nl and the diffusion coefficient, 
D^{t) = (1/2) dAnl/dt. The results for the latter read at short t 






2^-1 eVi 



(6) 



4 0F' 



and for large t 



D^{t)^^{A-n 



0r (1 - J)^/n 1 



(7) 
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The limiting value of the diffusion coefficient, Db(oo) ~ {v\y/eY is large com- 
pared to the diffusion coefficient of the one- dimensional random walk along 
the filament, i^b(O) = (1 — 7)/2. This broadening of the distribution occurs 
since the unbound motors lag behind the bound ones, which implies that also 
the rebinding motors lag behind those that have been bound for some time. 

The density profile on the filament can be evaluated analytically for large n 
and t. 



The exponential decay for n | v\,t expresses the large probability of unbinding 
for t ^ 1/e. Comparison with the results of Monte Carlo simulations shows 
very good agreement for times larger than about 8000 time steps [6] . 



3.2 Properties of the unbound motors 

Eventually every motor will unbind and diffuse in the surrounding fluid. We 
now discuss the effects of the filament on the behavior of the unbound motors. 

Position and longitudinal diffusion. One finds for the average velocity of 
unbound motors at small times -Uub ~ |'i^b(l ~ | ev^Trt) and at large times 
'i^ub(^) ~ '^h/W^^ \/t\- Whereas each individual motor has zero average velocity 
in the fluid, the statistical velocity v-ah is non-zero, since it is driven by unbind- 
ing from the cloud of motors moving on the filament. The longitudinal diffusion 
coefficient behaves at large times as D\\ = (tt — 2)v^/[7re^] -|- 1/4. The order of 
magnitude D\\ ~ "^b/^^ ieWs us that, like on the line, longitudinal diffusion is 
strongly enhanced by the unbinding from and binding back to the line. The 
perpendicular diffusion coefficient is modified much less, D±{t) ~ e v^/(2y^) 
at short times, while at large times D±{t) 1/4— l/(4e-\/7rt)- So the transverse 
diffusion starts out very small, and finally reaches its free space value. 

The spatio-temporal density profile of the unbound motors can also be derived. 



The presence of multiple tracks on microtubules can be modeled by internal 
states, which keeps the problem solvable [5,6]. 
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The same setup can be studied ind — 3 [5,6]. The Fourier-Laplace transformed 
density on the filament now reads 

" e + [3(1 -e)s + a (e"^- - 1) - b{e'- - 1)] I{r, s) " 



with a — |(e — 35) and b — |(6 — 67 — 35 — 5e) and involving the complete 
elliptic integral I(r, s) — (l/27r)^ Jq^ dqidq2 [3 + 3s — cosr — cosgfi — cosg2]~''^- 

The main difference is now that the density of motors on the filament decays 
faster, as ~ 1/t, because of a reduced return probability. Lack of space prevents 
us to discuss this any further, but we will present some results from continuum 
equations for the three-dimensional case below. 



4 Motor movements in confined geometries 



4-1 Scaling arguments and numerical results 



We have also studied these random walks in compartments with simple ge- 
ometries which are accessible to in vitro experiments [4]. In these geometries, 
a filament is immobilized to a surface, and the diffusion of unbound motors 
is restricted by confining walls. In the simplest case, the unbound motors can 
diffuse freely in the half space above the surface to which the filament is im- 
mobilized. By placing the filament in a quasi two-dimensional slab or in a 
rectangular or cylindrical tube, diffusion can be restricted along one or two 
dimensions perpendicular to the filament. We denote the linear extensions of 
the compartments by L±. Scaling arguments and simulations show that, at 
large times with i ^ i** ~ L\/ D^h, when the motors experience the pres- 
ence of the confining walls, these systems exhibit the same scaling behavior as 
systems with the same dimensionality, but without confining walls. For com- 
partments with unconfined diffusion in d± dimensions perpendicular to the 
filament, the effective velocity is given by v^^Nq with No{t) ~ t"'^-^/^. 



4-2 Continuum equations 



The continuum equations for the random walks of the motors can be solved in 
a similar way to the master equations of the lattice models. Since boundary 
conditions are easily implemented in the continuum model, we use them to 
derive analytical solutions for the half space, slab and tube geometry, see also 
Ref. [10]. We denote the coordinate parallel to the filament by x and the 
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coordinates perpendicular to it by y = (yi,y2)- The continuum equations are 
given by 



^ = L'ub Ap + 5(y) [eP - n^po] (11) 
ot 

dP dP ^ d^P _ 

with the bound-state velocity and diffusion coefficient and D\^, the unbound 
diffusion coefficient Dub, e = 2e/3, yfad = 27rad/3, and po{x) = i'^p{x, y = 0). 

Full three-dimensional space. Let us start with the full three-dimensional 
space. We use the Fourier-Laplace transformed form of Eqs. (11) and (12), 



sp{r, q, s) = -D^]:,{r'^ + ql + ql)p{r, q, s) + eP(r, s) - n^^dPoir, s) (13) 
sP(r, s) = l+ ivbrP{r, s) - D]3r'^P{r, s) - eP{r, s) + TTadPoir, s), (14) 

with the momentum q = (gi, q'2) and with the same initial conditions as before. 
This leads to 

00 _ 00 

f dqidq2 , . eP - n^po f .~ q 

= J 27 2^^^^' ^' = 2nD^, J "^'TTf- 

—00 



To obtain a convergent integral, we introduce a cutoff qc as the upper integra- 
tion limit. This cutoff corresponds to a certain range of attraction or thickness 
of the filament or a size of the motor and assures that the motor indeed returns 
to the filament. Defining 



^-/<l4T^ = i'"(^ + '«^' ^. = ^j^5. (16) 



we obtain pq = eP/ {2ttI ^ + Tfad) and 



P(r, s) = ( s - ivkV + D^r'^ e - J^^/^ _ ) 



V TT^\n{D^bqp ^) ) 

where the second expression is valid for small s and r. As before, the binding 
probability is obtained for r = 0. For sufficiently small s, the cutoff term can 
be neglected, and inverting the Laplace transform leads to iVo(t) = P{r ~ 
0,i) ~ 7fad/(47reDub^) in agreement with the result from the lattice model. 
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Half space. The case of the half space can be treated by the method of reflec- 
tion and superposition, sec e.g. Rcf. [11]. The part of the three-dimensional 
solution for 1/2 < can be reflected at the surface, so that the superposition 
^(2/2) 2/2) is obtained, which fulfills the boundary condition dp/ dy^ = at 
^2 = 0. The density profile is symmetric under refiection at the plane ^2 = 0, 
so we obtain just twice the solution for the full three-dimensional space and 
thus, also the double for the velocity and the displacement. 

Slab. Let us now consider a slab of height 2Lj_, where the filament is located 
at 2/1 = 1/2 = in the middle of the slab. The case where the filament is immo- 
bilized at the lower surface of the slab follows then by reflection and superpo- 
sition. The boundary conditions are dpjdy^ = at ^2 = i-^i. or, cquivalently, 
periodic boundary conditions at y-2 = ±L±. We can therefore expand the un- 
bound density into a Fourier series, p(. . . , ^2) = P(- • • > i.^j)e'^'^^^^ /{2L±) 
with Uj — Here and in the following we use the symbol u for discrete and 
q for continuous momentum variables. For the other variables, wc take the 
Fourier and Laplace transforms in the same way as above. Eq. (15) becomes 
now 

00 

Po(s,r)= j £ P{r,qi,ujj,s). (18) 



We integrate out qi as before. In addition, we approximate the summation 
over j by the term for j — and twice the integral 

S( ) = V 1 J do; 



where we introduced again a cutoff uic- We obtain 

4L_L-Dub/7rad 



P(r,s) 



-1 



s — ivy^r -|- e — e 1 + 



gub 



(20) 



for small r. For motors which have not yet reached the boundaries of the slab, 
i.e. for s :$> D.a^{T^ / Lj_y , we recover the result for the full three-dimensional 

space. For larger times or s <^ D^hiT^ / L^)"^ , the square root term is dominating 
in Eq. (20) and we obtain the two-dimensional result [10] with an effective 
reattachment rate 7fad/(2Lx). In summary, at large times the probability to 
be bound to the filament and, thus, the effective velocity i^b-^o behave as 



No{t) 



TTad/ {^T^^D^ht) for t<^U 

^tsd/{4:L±e^/^T^/D~^) for t :^ U 



(21) 
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with the crossover time i** = L^/ {nD-ah)- 

Tube. Finally, a tube with quadratic cross-section can be studied in the same 
way. In this case we have two discrete momenta uji^i and uj2,j corresponding to 
the transverse coordinates yi and y2, and Eq. (15) is now replaced by 

oo 

Po(g,r) = P(^>^M>^2j,s). (22) 

l^-^ij i,j=-oo 



As in the case of the slab, we approximate the double sum by the term for 
i — j — Q plus an integral with a cutoff ujc, 

^_ eP- TfadPo ( 1 , Ll J 27ra; do; ] 
- {2L^y D^^r^ + s + D^^n^ / D,b(r^ + c^) + J ' ^"""^ 



which leads to 

P(r,s) ^ {s -ivkr+ — ^ - + — ^In — ,^9 ) 24 



for small s and r. For s ^ Dub(7r/i^±)^, i.e., for motors which have not yet 
reached the boundaries, this expression yields again the solution for the three- 
dimensional case without confining walls, and for larger times, it leads to 



P{r,s) 



s{l + AL\—) -ivy,r 



n -1 



(25) 



and, therefore, to a constant probability to be bound to the filament, NQ[t) ~ 
1/(1 + 4L5^e/7rad), and a constant effective velocity v = fb/(l + 4L^e/7rad)- 
The crossover to the long-time regime is governed by the crossover time t** = 
V})/v X 7rad/(e47rDub) ~ -^i/(7r-Dub) where the last expression holds for large 
L± and agrees with the crossover time obtained for the slab. 



5 Summciry 



Molecular motors exhibit peculiar random walks which arise from the repeated 
binding to and unbinding from filaments. We have obtained analytical results 
for these random walks in systems without and with confining walls by solving 
the master equations of lattice models and the corresponding continuum equa- 
tions. Such random walks are relevant both for intracellular phenomena such 
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as, for example, the 'slow transport' in axons [12] and for artificial systems in 
nanotechnology in which molecular motors are used as transporters [13]. 

References 

[I] J. Howard, Mechanics of Motor Proteins and the Cytoskeleton, Sinauer 
Associates, Sunderland (Mass.), 2001. 

[2] M. Schliwa, G. Woehlke, Nature 422 (2003) 759-765. 

[3] A. Ajdari, Europhys. Lett. 31 (1995) 69-74. 

[4] R. Lipowsky, S. Klumpp, T. M. Nieuwenhuizen, Phys. Rev. Lett. 87 (2001) 
108101. 

[5] T. M. Nieuwenhuizen, S. Klumpp, R. Lipowsky, Europhys. Lett. 58 (2002) 468- 
474. 

[6] T. M. Nieuwenhuizen, S. Klumpp, R. Lipowsky, Phys. Rev. E 69 (2004) 061911. 
[7] O. V. Bychuk, B. O'Shaughnessy, Phys. Rev. Lett. 74 (1995) 1795-1798. 
[8] S. Klumpp, R. Lipowsky J. Stat. Phys. 113 (2003) 233-268. 
[9] S. Klumpp, R. Lipowsky, Europhys. Lett. 66 (2004) 90-96. 
[10] S. Klumpp, Ph.D. thesis, Universitat Potsdam (2003). 

[II] J. Crank, The Mathematics of Diffusion, 2nd Edition, Oxford University Press, 
Oxford, 1975. 

[12] L. Wang, C.-L. Ho, D. Sun, R. K. H. Liem, A. Brown, Nature Cell Biol. 2 (2000) 
137-141. 

[13] H. Hess, G. D. Bachand, V. Vogel, Chem. Eur. J. 10 (2004) 2110-2116. 



10 



